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SUMMARY 


The purpose of this study was to determine transient heat flow through 
ten exterior building walls for an air-conditioned building having reversible 
heat flow during a 24-hour period. The building walls were of various con- 
struction types; however, they each had essentially the same total heat 
transmission coefficient of about 0. 20 Btu/ft z -hr-° F. The inside room 
air temperature was maintained at a constant 70F. The variables involved 
in this study were the thermal capacitance of the walls and the distribu- 
tion of capacitance and resistance within the wall sections. 

A hydraulic analog computer was programmed to record transient 
temperature within the walls, and cumulative heat flow and instantaneous 
heat flux into the room. Three different daily cyclic input temperatures 
at the outside wall surface, derived from data for September at Wash- 
ington, D. C. , contributed one of the boundary conditions. A photographic 
technique described herein proved useful in recording results of the 
computer study. The ability of each wall to dampen load fluctuations 
was determined as was the time lag of heat flow into the room. A com- 
parison of the cumulative heat flow into the room during office hours 
was made for each type of exterior wall construction. This report con- 
tains: a description of the scaling technique used to program the com- 

puter (taking into consideration the time required for cancellation of 
transient ^ffects), the thermal properties of the ten wall sections, an 
approximate analytical technique for checking the accuracy of computer 
results, and results of programming the hydraulic analog computer 
for both a step and a sinusoidal input function for which analytical results 
are available. 

Results obtained in this study are given on figures and tables contained 
herein. Walls having high thermal capacitance result in a large thermal 
time lag, thus permitting the wall to store heat during the day and thereby 
delay full capacity operation of air conditioners until after working hours. 

A "peak load factor" is defined which provides a reasonably accurate 
relative measure of a wall’s ability to dampen heat flow variations in walls; 
it is tentatively concluded that a wall's damping characteristics are im- 
proved when the thermal resistance and capacitance are evenly distributed 
throughout the wall. The results of the study also showed that the greater 
the thermal capacitance, the greater the damping ability of the wall; 
this would be expected from theoretical considerations. 

It is concluded that the hydraulic analog computer can be successfully 
applied to solving heat flow problems subject to a few modifications 
necessary to reduce the average error of results. The influence of a 
building's location, the number of stories, required window space, and 
financial considerations governing construction are also factors in- 
volved in the selection of building walls in addition to air-conditioning 
requirements. It is recommended that the effects of direct radiation 
through windows, orientations of walls other than northerly, heat flow 
through the roof, and heat sources within the room be investigated more 
thoroughly before attempting to use the results of this study. 
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DEFINITION OF SYMBOLS 

Measurement 

Symbol 

Definition 

unit 

A 

Cross-sectional area of standpipes 

in. 2 

C 

Thermal capacitance 

Btu/ft 3 -°F 

c 

Thermal capacitance for a specified 
thickness of material 

Btu/ft 2 -° F 

C P 

Specific heat of solid 

Btu/lb-°F 

h 

Height of fluid in standpipes 

in. 

f 

Frequency of input 'temperature' 

Cycles /hr 

k 

Thermal conductivity 

Btu/°F-hr-ft 

P 

Density of solid 

lb /ft 3 

R 

a 

Analog resistance 

min/in. 2 

R 

P 

Prototype resistance 

°F-hr-ft/Btu 

t 

a 

Analog time 

min 

t 

P 

Prototype time 

hr 

T 

s 

Time scale = t /t 
a p 

min/hr 

U 

s 

Thermal energy scale = x V s J 

in. 3 /Btu ' 

V 

s 

Temperature scale 

in. /°F 

V 

Temperature 

° F 

4> 

Phase angle of input 'temperature* 

hr 

a, a 

Thermal diffusivity of solid 

ft 2 /hr 

Q 

Heat flow from inside surface to 

Btu/hr 


room air 


Q 

Average Q for one cycle 

Btu/hr 

P. F. 

Peak load factor = ^ 

unitless 


6 


X 

Distance through material 

ft 


HYDRAULIC ANALOG STUDY OF PERIODIC HEAT FLOW IN BUILDING WALLS 

by 

Rhode rick Hawk and William Lamb 
INTRODUCTION AND SCOPE 

As building technology progresses and comfortable working conditions become 
more desirable, accurate calculations of heat flow through building walls become in- 
creasingly important in selecting the proper air-conditioning capacity, particularly 
in office buildings. Determinations of heat gain allow one to choose the most suitable 
wall, both from economical and thermal considerations. 

This investigation was conducted to determine the periodic heat flow through ten 
typical building wall constructions (Table I) for an air-conditioned building. The 
problem was approached by using a hydraulic analogy method. However, these results 
do not represent the full picture of the air-conditioning problem since such important 
factors as direct radiation through windows,, infiltration into the room, heat sources 
within the room, heat flow through the roof, and orientations other than northerly were 
not considered. 

Although many investigations have been made in this field — by mathematical analy- 
sis, analog methods, and direct experimentation — it is felt that this study could be the 
beginning of a more advanced and complete study of the air-conditioning problem in 
office buildings. 

The model method used in this investigation is based on an analogy between heat 
conduction in materials and a fluid flowing in a network of small orifices and manometer- 
like tubes. The analogy between the hydraulic and heat phenomena is shown in the follow- 
ing table:* 


Temperature (°F) 

Thermal energy (Btu) 

Heat flow (Btu/hr) 

Thermal capacitance (Btu/°F) 
Thermal conductance (Btu/°F-hr) 


Height of fluid in a standpipe (in. ) 
Volume of fluid (in. 3 ) 

Fluid flow (in. 3 /min) 

Cross-sectional area of standpipe (in. 2 ) 

Volume of fluid passed through small 
orifice per unit time per inch difference 
in head of fluid in adjoining standpipes 
(in. 2 /min). 


POSTULATES OF THE TEST 

1) Three inputs were used for each wall. Two are typical of the 24-hour temper- 
ature ranges encountered in Washington, D. C. , during the air-conditioning month of 
September.! These are designated as inputs 1 and 2 and have minima at 60F and 75F, 
respectively. The third input is a temperature variation for one typical day during a 
heat wave, with a minimum at 82. 5F and a maximum at 100F. These inputs are shown 
in F igure 1 . 

2) All walls were designed for a total heat transmission coefficient (U) of approxi- 
mately 0.20 Btu/F-hr. These "U" values were made equal by the addition of selected 
insulation 1 in multiples of one-tenth of an inch. 

* For the derivation, see Appendix A. 

t September data were used for the inputs since data for the complete air-conditioning 
season (June, July, August, September) had hot arrived by the time the computer 
was ready for operation. 
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3) The room air was maintained at 70F. 

4) All walls were northerly exposures. 

5) The problems were solved for one-dimensional heat flow (no edge effects). 

6) The walls had no windows; there was no infiltration into the room; and there 

were no heat sources within the room. ( 

7) Moisture transfer within the walls was zero. 

8) . The only energy input to the external surface of the wall was by convection from 
the surrounding air. 

9) The outside air was varying in the periodic steady state. 

10) All resistances and capacitances were independent of time, temperature, or 
heat flow. 


Table I. Typical exterior building wall construction. 

U = . 20± 

Description 

4 in. CMU (light)*, Z in. cavity; 4 in. CMU (light); building 
paper 

4 in. brick; Z in. cavity; 4 in. CMU (heavy)f; insulation 

8 in. CMU (light); fin. furring; \ in. gypsum board; insulation 

4 in. brick; Z in. cavity; 8 in. CMU (light); insulation 

4 in. brick; Z in. cavity; 8 in. CMU (heavy); insulation 

"Sandwich" panel with Z in. glass fiber core clad in 18-gauge 
galvanized steel, to be installed in steel frame, primed and 
painted, inside and outside. 

8 -in. and 1Z in. concrete; f in. furring, \ in. gypsum board; 
insulation 

4 in. brick; 1 in. cavity; f in. plywood sheathing, Z x 4 studs, 

\ in. gypsum board; insulation- if needed 

Wood siding; f in. plywood sheathing, Z x 4 studs; -f- in. gyp- 
sum board; insulation 


* Hollow concrete block (shale aggregate). 

t Hollow concrete block (gravel aggregate). 

PROCEDURE 

Figure 2 shows how a thermal network with a finite number of segments was drawn 
directly from knowledge of the construction and the thermal properties of the materials 
used in the wall (Table II). The hydraulic analog network used in measuring the desired 
characteristics of the heat flow through the wall was constructed by appropriately scal- 
ing the problem. The procedures for scaling the hydraulic analog network, and the 
scaling sheets used for the actual calculations are given in Appendix A. 

The scaled hydraulic resistances* are set on the computer, and the appropriate 
standpipes (capacitors) (Fig. 3) are inserted in the places shown in the hydraulic 
analog network. 

* The hydraulic resistances were calibrated by the falling head method, described in 
reference 1. 


W all no. 
1 

Z 

3 

4 

5 

6 

7 & 8 
9 

10 
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STANDARD TIME (HOURS) 

Figure 1. 24 -hour temperatures for Washington, D. C. during September. 

The input function of time, f(t), is drawn on a sturdy paper to the appropriate 
scale, cut with scissors and made into a belt, and placed on a kymograph which feeds 
it at a scaled speed. The water level in a reservoir connected to the computer is 
made to follow this moving paper cam (by a system composed of contact switches and 
a motor), thus providing the input "temperature" (Fig. 4). 

The hydraulic apparatus representing the air-conditioner (and heater) is provided 
by a constant head reservoir placed at the end of the computer (Fig. 5). When water 
is flowing from the computer into the reservoir, it overflows into a calibrated tube. 

The rate at which the fluid flows is calculated by recording the height of the fluid in 
the calibrated tube at times representing hours in the thermal network. When the 
water is flowing from the reservoir into the computer, water is added to the reservoir 
by a system involving another calibrated tube filled with water, a valve, and a float and 
contact arrangement. Once again, the rate of flow is calculated by recording the height 
of fluid in this second calibrated tube. 

\ 

Although the data necessary for calculating heat gains were recorded manually 
while the wall sections were being run on the computer, it was decided that a method 
should be devised to record the temperatures in the various lumps of the walls (height 
of fluid in the capacitor tubes). A suitable method was obtained by using a photographic 
technique (App. F). 

After the computer components had been calibrated and the computer made ready 
for operation, two test problems with known solutions were run to test the operation 
and accuracy of the computer. The results of these two tests were reported previously 2 
and a brief resume is presented in Appendix D. 


* Numbers refer to references. 
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Table II. Thermal properties of the materials used in the typical wall constructions. 

\ 



Thermal 

Thermal 

Thermal 

Thermal 


conductivity 

capacitance 

diffusivity 

"wave 

length" 

Material 1 

K = 1 /R 

C - p x Cp 

a = k/C 

L = 0. 44 Ja/f * 

[Btu/°F-hr-ft). 

(Btu/ft 3 -°F) 

(Ft^/hr) 

m 

_ iisJ. 

Light CMUf 

0. 33 

14. 5 

0. 0228 

0. 325 

3.90 

Light CMU 4- air 

0.248 

3.6 

0. 0689 

0. 565 

6. 78 

Brick 

0. 75 

24. 6 

0. 0305 

0. 377 

4. 52 

Heavy CMUf 

1 . 00 

22.45 

0. 0445 

0.455 

5. 46 

Heavy CMU + air 

0.366. 

5.61 

0. 0653 

0. 550 

6.60 

F ur ring + air 

0. 06 

2. 14 

0. 0280 

0. 361 

4. 33 

F urring + fiber 
glass insulation 

0. 02Z5 

2. 4 

0. 0094 

0. 208 

2. 50 

Gypsum board 

0. 0925 

15.55 

0. 00595 

0. 166 

1. 99 

Fiber glass 
insulation 

0.021 

0.3 0 

0.-0700 

0. 580 

6. 95 

Steel sheet 

312.00 

58. 90 

5. 3 

4. 95 

59. 5 

Steel sheet + 
fiber glass 

1. 415 

0.787 

1. 114 

2. 27 

27. 2 

Concrete 

1.00 

22.45 

. 0. 0445 

0.471 

5. 65 

Plywood sheathing 

0. 067 

21. 5 

0. 00312 

0. 116 

1. 39 

2x4 studs + air 

0. 304 

2. 155 

0. 141 

0.810 

9. 73 

2 x 4 studs + fiber 
glass 

0.0225 

2.4 

0. 0094 

0. 208 

2. 50 

Wood siding 

0. 049 

17.4 

0. 00282 

0. 1145 

1.3 7 


* See Appendix B 
t Defined in Table I. 


HYDRAULIC ANALOG STUDY. OF PERIODIC HEAT FLOW IN BUILDING WALLS 


GLASS FIBER 

INSULATION / HEAVY CMU , AIR FILM 


OUTSIDE AIR 
(7.5 MPH WIND) 
TEMP - f (t) 




1 

vi 

r— 

AIR 



HEAVY 




l*. 





INSIDE AIR 
CONSTANT 
TEMP 


WALL SECTION 



16 


•r- 

i" 

2" 

16 

2 f 

9" 1" 

16 1 16 | 








■ 

D °F HR 
•' BTU 

0.29 

0.44 

0 86 

2.06 

0.07 

0.57 

0.07 

0.68 

7s btu 
Cp °F 

0 

7.50 

0 

0 

1.05 

1.17 

1.05 

0 


R 0, 


1 II! 

THERMAL PROPERTIES OF WALL 

Rc 

■JVr 

TCc 

. 

R °4 

Mn 

r b 

Rb ! Rb 
- w-pw - 

11 

nT 

: n I 

, i . 

Rb 1 Rq 2 Ri 

Rc 

it 

^( 0 * 03 ) 
—Mr 

"|Cc "pc 

@ 


1 

c B 

. 

L_ 

Cb 

1 

i. 




1 




i 1 

1 




THERMAL NETWORK TAKEN DIRECTLY FROM 
THERMAL PROPERTIES OF THE WALL. 


R| ^2 R3 ^4 ^5 ^6 R7 ^8 ^9 


1 ■ W — 
@) : 

-Mr 

-Mr 

-Mr 

— w 

rMr-] 

— Mn 

— H 


c T 

c r 

c r 

c 4 

C T" 


C7 

Ce 


HYDRAULIC ANALOG NETWORK DERIVED 
FROM THERMAL NETWORK 


Figure 2. Thermal properties arid analog circuits 
(example using wall 2). 








) 


HYDRAULIC ANALOG STUDY OF PERIODIC HEAT FLOW IN BUILDING WALLS 7 

RESULTS 

The cumulative heat flowing into the room as a function of time for one cycie for 
all ten walls and three inputs is plotted in Figures 6-27. The total heat gains for one 
cycle and for the working hours (7:00 a.m. -3:30 p.m. EST, 8:00 a. m. -4:30 p.m. EDT) 
are given in Table III. The cumulative heat curves for input 1 are given for only two 
walls since the heat flow into or out of the room is so small that accurate readings 
could not be obtained from the computer. This situation could have been remedied 
by rescaling the problems for input 1, but this would have added greatly to the length 
of work while increasing its value only slightly. 

The instantaneous heat flow to the room as a function of time for one cycle for all 
ten walls and three inputs is plotted in Figures 28-48. These curves were derived di- 
rectly from the cumulative heat curves. This method was considered best because 
the water (representing. Btu) did not always flow smoothly from the constant head reser- 
voir (air-conditioner). The time lags of heat flow into the room for the various walls 
with inputs 2 and 3 are recorded in Table III. The upper values in the table indicate 
the time lags from when the outside air temperature is at a maximum to when the in- 
stantaneous heat flow into the room is at a maximum. The lower values are the time 
lags for the minima. 


Table III. Heat gains and time lags for typical walls. 



Heat 

gain for 

Heat gain 

for office 

Time lag of heat 

Wall 

no. 

1 cycle (Btu/ft 2 ) 
input 2 input 3 

hours * 
input 2 

(Btu/ft 2 ) 
input 3 

flow to room 
input 2 input 3 

(hr)T 

Analytical 






max 

max 


1 

56. 8 

98. 2 

16. 8 

31.6 

4.42 

min 

3.75 

min 

* 3.52 






3. 00 

2. 00 


2 

56. 5 

99.8 

15. 1 

27. 1 

5. 00 
. 3. 00 

4.42 
4. 00 

5. 56 

3 

53. 5 

93. 5 

15. 7 

29. 6 

4. 00 
2. 00 

5. 42 
3. 00 

5. 02 

4 

49. 0 

86. 9 

14. 9 

29. 1 

9. 00 
5. 00 

7. 75 
7. 00 

5. 50 

5 

57. 4 

107. 0 

14. 0 

25. 9 

4. 00 
3. 00 

6. 42 
4. 67 

6. 61 

6 

54. 3 

92. 6 

24. 1 

40. 2 

1. 00 
1 . 00 

0. 00 
0. 30 

0. 52 

7 

53. 4 

93. 0 

14. 8 

2.8. 6 

7. 00 
3. 00 

6. 75 
2. 75 

5.10 

8 

61. 5 

108. 3 

16. 9 

30. 7 

9. 00 
7. 00 

7. 50 
6. 17 

5. 81 

9 

52. i 

89. 6 

14. 3 

27. 7 

5. 58 
3. 33 

4. 75 
4, 00 

5. 09 

10 

54. 0 

93. 6 

17. 7 

32. 0 

3. 75 
1. 00 

1. 25 
2. 50 

2.37 

* 7:00-3 

30 p.m. 

EST , 8: 

00-4:30 p. m 

. EDT. 





t Upper values indicate time lag between maximum outside air temperatures and 
maximum heat flow into room; lower values are time lags between minimum 
outside temperature and minimum heat flow to room. 
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In Figures 49 and 50, the cumulative heat curves and the instantaneous heat flow 
curves are given for wall 6, input 3, using three different mesh spacings of the therm- 
al network. A method of determining the necessary mesh size of a certain material to 
achieve 5% accuracy is described in Appendix B. The mesh sizes of the various wall 
materials determined by this method are given in the last column of Table II. 

In Figures 51, 52, and 53, the temperatures at various distances into wall 6 are 
given for one cycle for inputs 1, 2, and 3, respectively. In Figure 54, the energy 
stored in wall 6 throughout one cycle is plotted for all three inputs. These curves were 
drawn on the assumption that when the wall is at a uniform temperature of 70F, the 
energy stored in it is zero. This last group of Figures (51-54) is presented as a sam- 
ple of the information available on film for all ten walls and three inputs. 

In order to obtain an approximate check on the results of this study, the walls were 
divided into two segments, a sine wave was used as the input, and analytical solutions 
were obtained for the peak load factors and time lags ( the procedure used for the analy- 
tical solutions is described in Appendix E). The total resistance, capacitance, and 
density, and the simplified two-segment networks used in the approximation are 
recorded in Table IV for all of the walls. 

The concept of the peak load factor used in this report is a measure of the wall's 
ability to dampen load fluctuations and is defined in the following manner: 



The peak load factor is a ratio of the peak load to the average cooling load on the 
ai r- conditioner: 


P. F. = Q , /Q 

max avg 


The peak load factors obtained by both the analytical and computer solutions are 
recorded in Table V. The time lags obtained by both methods are recorded in Table 
IV. 
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Table IV. Thermal properties of the walls, 
and the two-segment simplified networks. 


WALL 

NUMBER 

TOTAL P 
POUNDS* 

TOTAL R 
F, HR/BTU* 

TOTAL C 

BTU/F* 

"TWO-SEGMENT" 

NETWORK 

l 

30. 85 

4.91 

4.22 


0.98 
2.11 = 

2.48 

= 

m 


2 

56. 06 

5. 00 

10. 78 


0.49 1 2*49 

7.5! T = 

2*02 
= 3.27 


3 

35. 03 

4. 88 

5. 49 

/ ' 

• . / 



I 

mZM 


4 

66.12 

5.03 

12. 20 


0.47 

7.50= 

2.88 

= u 

1.68 
= 4.70 


5 

80.31 

4.92 

14.77 


s 

1 1 

■M 


6 

5. 02 

5. 18 

0 . 62 


0.25 4.25 

031 X = 

0.68 
= 0.31 


7 

98.77 

5. 02 

15.80 


0.57 

15.0= 

2.05 

= 

2*40 
= 0.80 


8 


4.86 

23.30 


0.73 

22.5= 

r i 

2.17 
= 0.80 


9 

39. 00 

5.08 

9.93 



mm 

mmm 
m /v ,' 


10 

6. 29 

4.96 

3. 15 


wmjm 

Wmm 

mm 




* These values are for one square foot of the wall. 
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. DISCUSSION 


Total heat gains for one cycle 

For any given input, the total heat gains throughout one cycle should be the same 
for all ten walls. The proof of this is quite simple. The input function (temperature) 
is of the form A + B f(t), where f(t) is a periodic function. The thermal networks con- 
sidered are linear circuits, hence the heat flow to the room is of the form A/R + B/Z 
f(t), where _Z is the thermal impedance and R is the total resistance of the network. 

The integral of the heat flow to the room over one cycle will give the total heat gain. 

The integral of the periodic component over one cycle will always equal zero, and the 
integral of the non-time -varying component will equal (A/R) x (period). As A and the 
period are constants for any given input, and the total resistance, by design, is the same 
for all of the walls, the total heat gains for a particular input should be the same for all 
of the wall sections. 

An inspection of columns 1 and 2 in Table IV reveals that the heat gains for the 
walls are not the same, but vary by as much as ± 10%. The total resistance did vary 
somewhat from wall to wall. (The reasons for this and the other deviations are ex- 
plained in a subsequent section-"Sources of Error.") Closer agreement is observed 
if the heat gains are multiplied by the actual resistance of the wall (Table V). 

Heat flow to the room 

As noted above, the heat flow to the room is of the form A/R + B/Z f(t). Although 
the integral of this function for a given input will be constant for all of the walls, the 
amplitude and position of the heat flow curves will vary from wall to wall since the 
periodic component of the heat flow function has an amplitude and phase angle which 
depend upon the thermal impedance of the wall. It is in these heat flow curves that the 
thermal capacity of the walls becomes apparent. For purposes of comparison, a peak 
load factor, P.F. , was defined (page 8). This peak load factor provides a reasonably 
accurate relative measure of the wall's ability to dampen or attenuate the heat flow 
variation, since the difference in total heat gains will affect the P.F. values by only a 
small amount. (The difference in total heat gains is caused by a small shift - up or . 
down — of the complete heat flow curve. This shift will change the P. F. values by.a 
negligible amount). 

The ratio P.F. , used in this report, or some like measure, would appear to be a 
useful guide in determining the air-conditioning capacity necessary for a particular 
building. If the air-conditioning capacity were designed to meet the average heat flow, 
there would be a time during the day (most likely during the working hours) when the 
air-conditioner would not be able to keep up with the heat flow into the room. And the 
larger the P.F. value for that particular wall, the more the temperature in the room 
would rise when the air-conditioner lagged behind the heat flow. It is obvious, then, 
that the smaller the P.F. value, the smaller the necessary air-conditioning capacity, 
and, consequently, the smaller the costs of installing and operating the air-conditioners. 
It is observed from Table V that there is not always close agreement between the 
analytical and computed P.F. values. This is to be expected since the input used for 
the analytical approximation was a sine wave, while the input curve .used on the com- 
puter is represented by a Fourier series. But what is important is the relative agree- 
ment between the various walls. Three of the best walls by both methods are walls 5, 7, 
and 8. These are also the walls' with the largest thermal capacities. This correlation 
between high damping power and large thermal capacitance is in good agreement with 
theoretical expectations since the damping characteristic depends upon the impedance, 
which in turn depends upon the capacitance. However, the P.F. value is only approxi- 
mately a function of the thermal capacitance and depends to a certain extent on how the 
resistance and capacitance are distributed throughout the wall. 

The dependence of the P.F. value on the distribution of the resistance and capacit- 
ance is demonstrated by the results obtained for walls 5 and 7. Wall 7 has a large 
capacitance, but it is concentrated almost entirely at the outside of the wall. Although 
wall 5 has a slightly smaller capacitance, it is divided almost exactly in half with the 
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Product 

Table V. Relative heat gains and 
peak load factors for typical walls. 

of heat gain for one Peak load factor. 

,Q 

P-F. =(g=^ 

cycle and total R (°F-hr 
input Z A*(%) input 3 

-ft 2 ) 

A*(%) 

input Z 

input 3 

VU avg ' 
Analytical 

Z8Z.0 

+4. 1 

488. 0 

+ Z. 1 

1 . 49 

1 . 36 

1.25 

Z8Z.0 

+ 4. 1 

499. 0 

+4. 3 

1.34 

1.Z7 

1.16 

Z61.0 

-3. 7 

457. 0 

-4. 3 

1.35 

1.Z8 

1.18 

Z55.0 

-5. 9 

450. 0 

-5. 9 

1. 18 

1.16 

1.17 

Z71.0 

0. 0 

5Z6. 0 

+ 10. 

1.16 

1. 19 

1.12 

Z8Z.0 

+ 4. 1 

480. 0 

+ 0. 4 

1.66 

1. 47 

1.. 31 

Z68.0 

-1. 1 

467. 0 

-Z. 3 

1 . Z5 

1.Z4 

1 . 14 

Z74.0 

+ 1. 1 

48Z.0 

+ 0. 8 

1.17 

1. 1Z 

1. 08 

Z66.0 

-1.9 

456. 0 ( 

-4. 6 

1.40 

1.26 

1. ZO 

Z68.0 

-1.1 

464. 0 

-Z. 9 

1 . 4Z 

1.41 

1.Z7 


Wall 
no. 

1 

Z 

3 

4 

5 

6 

7 

8 

9 

10 

= Deviation from average 


(Average = Z71 for input Z 
= 478 for input 3) 


bulk of the resistance coming in the central part of the wall (Table III). The P. F. 
values are lower for wall 5 by both the computer results and the analytical approxima - 
tion. The conclusion then is that a wall's damping characteristics are improved when 
its resistance and capacitance are distributed more evenly throughout the wall. But 
this conclusion is only approximate, and the distribution of resistance and capacitance 
should be more thoroughly investigated because of its importance in air-conditioning 
design. 

Time lag of heat flow 

The time lag is larger at the maximum than at the minimum in almost every case 
(Table IV). This difference is easily explained by realizing that the time lag is a 
function of the impedance which is partly dependent upon the frequency of the input. 
Since the input is not a true sine wave, the frequency of the periodic component varies 
between the maximum and the minimum and is largest when the input curve is at its 
maximum (this is borne out by a close inspection of the input curves - Fig. l). Thus 
the impedance — and therefore the time lag — is greatest at this point. 

If the time lag is large enough that the peak heat flow occurs well after working 
hours, some of the stored heat will flow back out into the cooling night air instead of 
into the room. This has an important consequence on the air-conditioning capacity 
in either one of two ways. 

The first case is where the air-conditioners only operate during the working hours 
and it is assumed that the room and the walls will cool down naturally during the night. 
(This principle was first used by the Spanish settlers in Southwestern United States — 
Albuquerque, N.M. , for example — who built adobe houses with walls Z or 3 feet thick. 
During the cooler desert nights, they ventilated their homes and cooled the walls from 
both sides.) In this situation the air-conditioners need only be designed to meet the 
largest demand during the working hours, which will be in the vicinity of the average 
heat flow for the cycle, and possibly even lower. This fact has great .economic sig- 
nificance, both in the initial investment and in the costs of operation and maintenance. 
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The second case occurs most frequently in areas where the seasonal temperatures 
are unusually high. Here the walls and the room are not able to cool naturally during 
the night, and it is necessary to keep the air-conditioners running constantly through- 
out the cycle. Although it is necessary to design the air-conditioners to meet the peak 
load, they will be running at their full capacity only outside of working hours. This 
consideration is especially important in large industrial areas where the power con- 
sumption during the day is already at a point near the upper limit of the power supply. 

Table IV shows the analytical and computed values of the time lags. Heat gains 
for office hours are recorded in the same table. (These values have been adjusted 
assuming all of the heat gain'* for one cycle to be equal.) It is observed that walls 
7 and 8 have large time lags, ">ut for heat gains during the working hours, wall 7 is 
one of the best of all the walls, while wall 8 is among the worst. This is the result 
of a cancellation effect between the peak load factor and the time lag. When the time 
lag is large — on the order of 6 hours or more — the average heat flow to the room 
during working hours is less than the average for the cycle. (For a large time lag, the 
heat flow to the room is near the minimum in its cycle during the working hours. See 
Figure 4, for example.) If the heat flow to the room during the working hours is less 
than the average for the cycle, it becomes obvious that the more the heat flow is 
attenuated the greater the heat gain will be during the office hours. This was the case 
for wall 8, as can be seen by comparing its relative heat gain and peak load factor 
with those of wall 7. 

It might be pointed out that wall 6 has the smallest time lag and the greatest peak 
load factor (thereby the least attenuated heat flow) and the greatest heat gain for the 
office hours. This is a direct consequence of this wall's extremely low thermal capaci- 
tance. 

Sources of error 

As noted earlier, the errors in total heat gains run as high as ±6-7%. Errors of 
this magnitude were not expected in view of the results of the test problems (App.D). 

To account for these errors, an appraisal was made of the maximum error contributed 
by each of the components in the computer. 

The hydraulic resistance is directly proportional to the viscosity of the fluid, and 
the viscosity is dependent upon the fluid temperature 3 . At the time the tests were run, 
the air temperature in the laboratory varied throughout the day by as much as ±6C. 
Variations of this magnitude in the fluid temperature during a test would result in 
errors of 1 0% or more. But every effort was made during the tests to keep the resistors 
set within ±1C of the fluid temperature. Variations of 1C contribute an error of ±2% 
in the overall resistance of the wall. 

The period of the input curve varied by about ±0. 6% from that used to . scale the 
problems. On wall 5 this error was about ±1%. It is believed that this error is caused 
partly by the kymograph's speed failing to remain constant, but more by a slight slippage 
of the cam (input curve). This error in the period affects the resistors directly and 
changes the time lags by about ±15 minutes. 

Both the input reservoir (outside air temperature) and the constant head reservoir 
(air-conditioner and heater) were set on the computer with respect to the reference lines. 
An error of approximately ±0. 6% is possible from the setting of each of these two com- 
ponents. 

An error of roughly ±1% is a reasonable estimate for the calibration and setting 
of the hydraulic resistors. Also included in this estimate would be the effect of any 
unaccounted-for resistance in the stopcock valves. 

An error of ±1% results from the calibration and selection of the correct capacitor 
tubes. 

As the programmer followed the input curve in small steps up or down, there was 
a certain amount of induced error. This possibility that a given input curve is not 
followed in exactly the same manner in two tests accounts for an error of ±1%. 
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The method used to determine when the effect of transients was insignificant (App. 
C) should hot allow for an error greater than ±0. 5%. 

The reading of the calibrated tubes for the constant head reservoir, and the method 
used to record the data from the films (Fig. 51-54) could produce ari error of about 
± 1 %. 

These sources of error account for a maximum possible error of approximately 
±8%. As can be seen in Table V only a few of the deviations approach this magnitude, 
and most of the errors are much smaller. This is to be expected since a certain 
amount of cancellation of errors will occur during the normal operation of the com- 
puter. The method generally used to calculate the average error consists of taking 
the square root of the sum of the squares of the maximum possible individual errors. 

By this method, the predicted average computer error is 3. 1%. The average of the 
computer errors given in Table V is 3.2%. 

It was concluded from this agreement between the predicted and the actual errors 
that the hydraulic model method can be successfully applied to solving heat transfer 
problems, but that every effort should be made to reduce the average error. Many 
of the inaccuracies causing this error could be eliminated with a few minor improve- 
ments to the computer. Immersion of the computer in a constant temperature bath 
would eliminate the largest source of error. The design of new resistors with greater 
resistance ranges would eliminate the necessity of occasionally using several resistors 
in series, which causes a certain unknown amount of error. 

. . ; I 

CONCLUSIONS AND RECOMMENDATIONS 

Before the results of this investigation, or any other study like it, can be used to 
select the proper building wall or air-conditioner, it must be remembered that many 
other factors influence their selection; for instance, the location of the building site, 
number of stories, window space desired, the financial conditions governing its con- 
struction, and the personal preferences of the prospective owner. 

We have seen that walls 7 and 8 have very large thermal capacities which gives 
them low peak load factors and large time lags, desirable characteristics for the 
walls of an air-conditioned building. But these characteristics are the result of the 
fact that the exterior of these two walls consists of 8 in. and 12 in. poured concrete, 
respectively. It is not only obvious that walls constructed in this manner will be very 
expensive to build, but that they will also be very impracticable walls in a building 
much higher than two stories. On the other hand, if the building is to be two stories 
high and air-conditioning is one of the major problems, these walls might be the 
most suitable — both economically and from a thermal standpoint. These walls might 
also be satisfactory if they were to be used in the , first two stories of a multi-storied 
skyscraper. 'Heavy 1 walls such as walls 7 and 8 would provide the strength needed 
in the base walls of the skyscraper. 

Another important factor for consideration is the relative costs of the materials 
and construction of these walls. Considering complete costs of construction, walls 
7 and 8 are the most expensive, wall 1 0 is by far the least expensive, and the remainder 
of the walls fall in approximately the same price category. Wall 6 is very easy and 
inexpensive to erect, but the cost of the metal panel offsets this saving. Although 
wall 1 0 is a very inexpensive wall to build, it would very often be considered an unde- 
sirable wall for buildings other than residential structures. Reasons for this would 
include maintenance costs, general appearance, and the fact that it would not be a suit- 
able wall for multi-storied buildings. 

Relative ease in framing windows and maintenance costs are two other important 
factors for consideration. Where a great deal of window space is desired, it would 
be an easy and inexpensive operation for walls 6 and 10, but quite expensive for the 
others. 'Where maintenance costs are considered, it is found that wall 6 is more 
satisfactory than the others. It would also be a simple procedure to dismantle wall 6 

and reassemble it in another location - a characteristic not possessed by the other 
walls. 
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In considering relative costs and thermal characteristics of these walls, it should 
be remembered that this study was subject to many limitations and simplifying 
assumptions. In the first place, cost comparisons are difficult since many of these 
walls are unrealistic constructions. The overall transmission coefficient was made 
a constant value for the various walls by adding successive layers of insulation in 
thicknesses of one-tenth of an inch. In reality, insulation is seldom used in thick- 
nesses less than one inch. To cut it to a smaller thickness is a difficult and eco- 
nomically unfeasible process. 

Factors not included in this study include direct radiation through windows, 
orientations other than northerly, heat flow through roof, and heat sources within 
the room. Actual heat conduction through the walls makes up only a small fraction 
of the total heat gain considering these factors. In particular, direct radiation through 
windows will have a significant effect on the time lags and peak load factors achieved 
in this report. A more thorough investigation of these considerations should be made 
before any attempt is made to use the results of this study in selecting the most suit- 
able wall or air-conditioning capacity for a particular building. 
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Figure 11. Cumulative heat passed through wall 3, input 3. 
Total Btu/day-ft 2 ^93.5. 
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Figure 16. Cumulative heat passed through wall 6, input 1. 
Total Btu/day-ft 2 ^ -4. 9. 





HYDRAULIC ANALOG STUDY OF PERIODIC HEAT FLOW IN BUILDING WALLS 


3 

I- 

CQ 



Figure 18. Cumulative heat passed through wall 6, input 3. 
Total Btu/day-ft 2 ^ 92. 6. 



Figure 19. Cumulative heat passed through wall 7, input 2. 
Total Btu/day-ft 2 ^ 53.4. 
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Figure 23. Cumulative heat passed through wall 9, input 2. 
Total Btu/day-ft 2 ^ 52.4. 
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Figure 26. Cumulative heat passed through wall 10, input 2. 
Total Btu/day-ft 2 ^ 54. 0. 
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Figure 27. Cumulative heat passed through wall 10, input 3. 
Total Btu/day-ft 2 ^ 93.6. 
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Figure 31. Heat flow through wall 2, input 3. 
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Figure 35. Heat flow through wall 4, input 3. 



Figure 36. Heat flow through wall 5, input 2. 
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Figure 37. Heat flow through wall 5, input 3. 
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Figure 40. Heat flow through wall 6, input 3. 



Figure 41. Heat flow through wall 7, input 2. 
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Figure 42. Heat flow through wall 7, input 3. 






HEAT FLOW INTO ROOM HEAT FLOW INTO ROOM 

(BTU/HR.) (BTU/HR. ) 










34 






HEAT FLOW INTO ROOM (BTU/HR. ) 


HYDRAULIC ANALOG STUDY OF PERIODIC HEAT FLOW IN BUILDING WALLS 



SUN TIME (HOURS) 

Figure 49. Cumulative heat curves for wall 6 
with 3 mesh sizes 



Figure 50. Heat flow curves for wall 6 and 3 mesh sizes. 
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SCALING THERMAL NETWORK FOR THE 

\ 

HYDRAULIC ANALOG NETWORK 

If we consider heat transfer through any material in one direction only, the field 
equation is given by: 


8 z v _ 

8X 2 



k at 

P 


where, v = temperature difference, °F 
x = distance through material, ft 
p = density, lb/ ft 3 
( Cp = specific heat, Btu/lb°F 

k = thermal conductivity, Btu/°F-hr-ft 
tp = prototype time, hr 


A finite difference equation can be derived from the field equation as follows: 


(Al) 


v (x + Ax) - v (x) v (x) - v (x-Ax) _ 9v (x) 

(Ax) 2 /k (Ax) 2 /k P 9t p 

If we note that (Ax/k) is the thermal resistance (multiplied by unit area) of a specified 
thickness of material, and that Ax • pC p is the thermal capacitance (per unit area) of 
the same specified thickness of material, we have 

vi-vp _ v o - v -i - q 9 v o # (A2) 


Equation A2 is an approximate solution for v 0 in the thermal network shown in 
Figure Al. The finite difference equation used for the solution of h 0 on the hydraulic 
analog network (Fig. A2) is given by: 

hi - h 0 h 0 - h-i _ A 9h 0 (A3) 


where, A = area of standpipe, in. 2 

h = height of fluid in standpipe, in. 
R a = analog resistance, min/in. 2 
t a = analog time, min. 


It is immediately obvious that equations A2 and A3 are analogous, and that we are now 
ready to scale the thermal network for the hydraulic analog network. 

The first step is to divide the wall into a finite number of segments, taking into 
account the number of available positions on the computer. In general, the thickness 
of the segments should be smaller where the temperature gradient is large. (In this 
study, the temperature gradient was largest in the outer portion of the wall. ) Once 
the thicknesses of the segments have been determined, the thermal capacitance for each 
segment is determined (see Figure A3 for a sample problem): 


Thermal capacitance = Axp Cp - C> 


Btu 


(unit area). 


A 2 
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Figure Al. Thermal network. 


Figure A2. Hydraulic network. 


As seen from equations A2 and A3, the hydraulic analog of thermal capacitance 
is A, the cross-sectional area of the standpipes. A ratio of the standpipe area to the 
thermal capacitance is chosen: 


A 

C 


in. 


2 o 


Btu 


= constant. 


(A4) 


The value of this ratio is determined by the areas of the available standpipes, and by 
the fact that this ratio is a constant for all of the segments in any one wall problem. 

From the two finite difference equations, it is observed that the hydraulic analog 
of temperature difference is the hydraulic head. A temperature scale V s , which is a 
ratio of the hydraulic head to the temperature difference, is chosen taking into account 
the limits of the computer (i.e. , the temperature range of the input curve and the height 
of the standpipes): 

V = h/v, in. /°F, (A5) 


must 


A time scale T s , the ratio between the analog time and the prototype time, .. 
be selected, remembering the limits of the input curve and the kymograph speed, and 
checking to see that it does not go above or below the limits of the linear range of the 
hydraulic resistors: 


T = t /t , min/hr. 
s a p 


(A6) 


The thermal network resistances must be determined next. The procedure for 
making this calculation is demonstrated in the. sample problem in Figure A3. The 
analog resistances are then scaled from the thermal resistances according to the 
following equation: 


R a - T g (C/A) 11^, min/in. 2 . 


This last equation is derived by substituting the three scaling ratios already obtained 
eq A4, A5, and A6 — into the finite difference equation for the thermal network: 

Starting with = c .§£&. 


R 


at 


we multiply through by V 
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hi - hg hg - h_i — r~ 3^o 

R . R at 

P P P 


divide by T : 

. s 


hi - hg hg - h _ 1 ^ _ 9hg 

T R ‘ T R at 

s p s p a 


and multiply by A/C: 

h-i - hg hg - h_i _ A _ahg_ 

T s (C/ A )R p T s (C/ A )R p 

A comparison of this last equation with equation A3 reveals that, 

R = T (C/ a )R , min/in. 2 . . 
a s ' A p 

The ratios A4 and A5 determine U s , the energy scale, which is a ratio between 
the volume of water and the thermal energy: 

U = (A /-^t) V , in. 3 /Btu. 
s ' C s 

The input curve is then scaled in accordance with T s , V s and the speed of the kymo- 
graph. The constant head reservoir (air-conditioner and heater) is set at the proper 
height on the computer with respect to the input curve, and the computer is ready for 
operation. 





STEP ONE 


STEP TWO 


STEP THREE 


Figure A3. Sample problem of resistance 
and capacitance scaling procedures. 


= 6 Btu/ ft 3 -° F. 


= Ax x /k x 
= Ax 2 /k x 
= Ax 3 /k x 


1 ft 2 = 

1 .ft 2 = 

1 ft 2 = 


3° 

F-hr /Btu 

Ci = Ax • 

1 ft 2 Cx = 3 °F/ Btu 

2° 

F-hr / Btu 

C_ 2 = Ax 2 

1 ft 2 C x = 2° F / Btu 

1 . 

5° F-hr/ Btu 

C 3 = Ax 3 

• 1 ft 2 C x = 1. 5° F/ Btu 
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HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


PROBLEM! _HEAI_FJ 

TIME SCALE: I MIN: L 

WALL PROFILE 


OUTSIDE SURFACE 
(7,5 MPH WIND) 

SECTION I LIGHT 
CMU 9/16" 


$ S ^ 

i 

JM 


SECTION 3 LIGHT 
CMU 9/16" 


AIR SPACE 
~ 2 " 


BUILDING PAPER 
1 2 LAYERS 

SECTION; I LIGHT 
CMU 9/16" 


I ^ 

£ > 

L> «Vl 


SECTION 3 LIGHT 
CMU 9/16" 

INSIDE SURFACE 
(STILL AIR) 


test no':* L 


I »F: 0.4 INCH 

I 


COMPUTED BV R- H, 

BTU: 0.0 436 INCH 3 I °F HR/RTII: 9.17 MtN/IN 2 


STANDPIPE 

AREA 

(INCH 2 ) 

LUMP 

THICKNESS 

(INCHES) 



0.0740 

0.563 

0.0407 

1.250 

0.041 1 

1.250 

0.0741 

0.563 

0 

2.000 

0.'04l 1 

0.563 

0.0408 

1.250 

0.0410 

1 .250 

0.04l 1 

0.563 


CENTERS 


BTU IN 2 (INCHES) 


0.425 3.89 8.50 

0.906 0.385 3.53 9.05 

1.250 0.420 3.85 8.00 

/ 

0.906 0.385 3.53 8.70 


0.906 0.385 3.53 13.15 

1.250 0.420 3.85 11.75 

0.906 0.385 3.53 9.05 

" 0.855 04 7784" 


* Test no. indicates wall no. 


ROOM AIR TEMPERATURE = 26.3 °C 
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HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


PRO 8 L E M : HEAT FLOW THROUGH WALLS 


TIME scale: IMIN = _J HR 


WALL PROFILE 


OUTSIDE SURFACE 
(7.5 MPH WIND) 


I °F - A— INCH 


4 FACE BRICK 
(3-5/fl") 


7.50 l 0.44 


AIR SPACE 

1.5" 


0 % 


BUILDING PAPERS) T* 

_J GLASS FIBER H— 

INSULATION 0 5" 0.013 2.06 

SjpiON ' HEAVY, | 053 q 07 

V. 

L, ^ * 

§ * S 1 1 .170 0.57 

5 S «Sl 

N - 


SECTION 3 HEAVY 

CMU S/16 

INSIDE SURFACE 
(STILL AIR) 


1.053 0.07 


TEST NO.; 2 

date: 21 MA RCH 1962 

COMPUTED BY R H. 

i BTU ; 0 02 75 inch 3 i°F HP/RTU: 14.55 min/in 2 
standpipe; lump ; INCHES fl r ~ 0 r; 

AREA j THICKNESS T0 o F hr ^ SETTING 

(INCH 2 ) j (INCHES) CENTERS BTU | N 2 (INCHES) 



— 

A OAC 

• 

0.305 

4.44 

7.23 

O.IZSc 

U.SUO 






A AAfi 

0.906 

0.11 

1.60 

12.4 

O.IZB7 

V.9VV 





A 90 OO 

i 

f\ QAC . 

0.906 

0.11 

1.60 

12.45 

U.leo? 

u.yuo 







0.906 

0.1 1 

1.60 

12.45 


0.1285 0.906 — 

‘0 1.50 2.734 j 3.010 

I 

0 0.5 | 

O 0722 0.563 j 


10 7.53 

10 4.72 


13.8 6 9 

10 4.9 


0.0400 1.25 

0.0400 1.25 

0.0723 0.563 


0.906 i 0.178 2.59 13.82 

1.25 0.285 4.15 11.3 

0.906 0.178 2.59 10.75 

0.715 10.40 11.6 


1=5.00 U=.200 


ROOM AIR TEMPERATURE = 27.5 °C 


$LASS FIBER INSULATION 0.2 
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HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


TIME SCALE: I MIN: 


WALL PROFILE 


OUTSIDE SURFACE 
(7.5 MPH WIND) 

SECTION I LIGHT 


°F: Q,4 INCH 


I RTU- 0.074 
[standpipe I LUMP 


DAT E I 21 MARCH 1 962 

COMPUTED BY' R.H. 

■I 3 l°F HR/RTU; ,5.72 MIN/IN* 


iTANDPIPE LUMP iwrHF*? R„ R R ! 

AREA THICKNESS to 0p > R ^ SETT ° |NG 

(INCH 2 ) (INCHES) CENTERS BTU IN 2 (INCHES): 


0.300 | 1 .313 


0.50 2.86 110.34 


0.0702 1. 25 


0.0708 1.25 


1.282 10.375 2.12 11.43 


1.25 i 0.25 1.43 12.77 


1.00 > 2.00 


1.25 0.25 1.43 12.76 


0.0710 1.25 


1.25 10.25 1.43 14.6 


SECTION 3 LIGHT 

CMU 1-5/16“ | 59 

PAPER COVERED -2 

AIR FURRING SpACE 0.098 


rad 3/4 "» 12 " h 

GYPSUM BOARO 
1 / 2 " 

INSIDE SURFACE 
; (STILL AIR ) 


0.040- 


0.0716 1.25 


0.297 1.313 


0 0261 0.75 

0.1184 0.50 


1.282 0.375 2.12 14.15 

1.032 1.005 5.75 9.45 

.625 0.970 5.55 4.95 

0.903 5.16 10.6 


ROOM AIR TEMPERATURE = 27. 5 °C 
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APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


ppnm fm- HEAT FLOW THROUGH WALL 


TIME scale: I MIN-__1__HR 


°F : 0.4 INCH 


C R 

WALL PROFILE. BTU °F HR 


OUTSIDE SURFACE 
(7.9 MPH WIND) , 


4" FACE. BRICK 
( 3-9/8“ ) 


7.50 0.44 


AIR SPACE 
1 . 8 " 


SECTION I LIGHT 
CMU 

( 1-5/16) 


Qj 

0.86 

0,005 0.80 - 

1.59 0.50 


& £ 
5 £ 


| * V * 


* S 

*" P 

o> ki 

** 


SECTION 3 LIGHT 
CMU 

(1-5/16) 

INSIDE SURFACE 
(STILL AIR) 


1.59 0.50 


TEST NO.' : — 

DAT E . 2 0 MARCH 1962 

compute^ ry: R.H. 

R T U : 0.02 58 INCH 3 i°F HR.'RTH: 15.90 MIN/IN 2 


(INCH 2 ) : (INCHES) 


* 2 r S 

I ^ | n I 162 100 / 2.00 


215*5.03 U**. 199 


- — 1 

0.1166 

— ! 

0.906 j 

i 

0.1170 

0.906 

> 

0.1169 

0.906 

0.1167 

0 906 

. 

0 

— o — 

• 

i 

! 

1.80 

- 0 . 20 — 

0.1001 

1.313 

0.0238 

1.25 

0.0244 

1.25 

0.0241 

1.25 

0.0240 

1.25 

0.1001 

1 .313 




— 1 

INCHES 
TO ' ! 

CENTERS j 

r p 

°F HR 

"btu 

j.\ R o. 

i MIN 
IN 2 

I settin 

(INCHE! 

| 

0.305 

6.14 

8.15 

10.95 

0.906 

i 

i 0.11 

2.21 

10.9 

1 

0.906 

i 

1 ■ 

0.11 

2 21 

11.18 

i 

I 

0.906 

1 

an 

1 

2.21 

10.50 

. i 

; 

i.’’ , 

12.00 

i 

5.85 

3.109 ! 

i 

1.965 

9.00 

5.95 

l 

1 

t 

• 1 
■ ■■ i 

l 

! 

; 

i - - - 

10.20 

10.87 

i 

1 .283 

! ' 

0.375 

5 96 

ll.l 

1.25 

0.2$ 

3.9T 

11.47 

1.25 

■i 

0.25 

i . 

3.97 

8.05 

1.25 

0.25 

3.97 

11.87 

1.283 

0.375 

5.96 

15.8 


0.93 

14.80 

5.00 


ROOM AIR TEMPERATURE = 27.5 °C 
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APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


prori fm: HEAT FLOW THROUGH WA1 


TIME SCALE: I MIN: I HR 


WALL PROFILE 


OUTSIDE SURFACE 
(7.5 MPH WIND) 


l°F: 0.4 INCH 

R I 


4 FACE BRICK 
(3-5/8") 


7.50 0.44 


1.9 

AIR SPACE 




% 

2.04 

O l“ GLASS FIBER j 
INSUL. W/POUSHED 1 


-0.0025 

ro.40 

ALUM COATED PAPER 1 
BOTH SIDES | 


2.46 

• 

0.18 

SECTION 1 HEAVY 
CMU 1-5/16" 


TEST NO •_ 5 

DAT E : 20 MARCH 1962 

COMPUTED bt: R.H. 


<M ^ 


SECTION 3 HEAVY 
CMU I -5/16" 

INSIDE SURFACE 
(STILL AIR ) 


2.46 

0 ~ 


IBTU: 

0.01785 

.INCH 3 

l°F HR/BTU i - 22.5 

MIN/IN 2 

STANDPIPE 

LUMP 

INCHES 

R P 

R 0 

i R ° 

AREA 

THICKNESS 

TO 

°F HR 

MIN 

SETTING 

(INCH 2 ) 

(INCHES) 

CENTERS 

BTU 

IN 2 

l(INCHES) 


0.839 0.906 : : 

0.906 0.11 2.48 11.13 

0.837 0.906 

0.906 0.11 2.48 10.95 

0.838 0.906 

0.906 0.1 1 2.48 10.95 

0.840 0.906 

4.85 

4.78 

0 I 90 3.11 2.585 58.20 4.8 

4.52 

— 0 0.10 - 

9.7 

0.1102 1.313 

~~ 1.54 0.222 5.00 12.14 

0.0368 1.76 — — 

1.47 0.220 4.95 10.8 


0.0245 1.17 


0.0247 1.18 


0.0183 0.88 


0.1107 1.313 


1.18 0.176 3.96 8.55 

r 

1.03 0.155 3.48 10.2 

MO 0.156 3.51 13.4 

.656 0.774 17.40 4.4 


IR=4.92 


U = .203 * 


ROOM AIR TEMPERATURE = 24.0 °C 



GLASS FIBER INSULATION 2 


APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


ppr>Ri pm- HEAT FLOW THROUGH WALLS 


TIME SCALE: I MIN 


WALL PROFILE 



) date: 

COMPUTED BY! 

|0 F; _04_INCH ' IBTU: -296 INCH 3 I °F HR/BTU : _3j»6__MIN/IN 2 

STANDPIPE j LUMP INCHES I ^ K ^0 

AREA j THICKNESS T0 op HR MIN SETTING 

(INCH 2 ) ! (INCHES) CENTERS BTU IN 2 (INCHES) 


25 4 


0.275 0.531 1.84 12.1 



0.1815 

0.05 

.0241 

0.50 

.0240 

0.50 

j 

.0244 

0.50 

.0245 

0.50 

.1812 

0.05 


0.50 I 1.061 3.67 8.52 


0.50 1.061 3.67 8.57 


0.50 1.061 3.67 12.42 




£R = 5.175 U = 0.1966 ; ROOM AIR TEMPERATURE = 25.5 °C 












































APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET- 


PROBI fm: HEAT FLOW THROUGH WALLS 


TIME scale: I MIM = I HR 


WALL PROFILE 


I °F r _0-4 INCH IBTU: O H* 


DATE . io AU GUS T 1962 

COMPUTEO RY- H.H. 

l°F HR/BTLJr _ 3.6 MIM/IM 2 



XR = 5.175 U =01966 


ROOM AIR TEMPERATURE = 24.0 °C 



























GLASS FIBER INSULATION 2 


APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


PRO B L E m: HEAT FLOW T HROUGH WAj 


' I ME scale: I MIN: !_t HR 


WALL PROFILE 

OUTSIDE SURFACE 
(7.5 MPH WIND) 
PAINTED-® 


i °e - 0X54. inch 


0 | 0.25 

0.245 j 0.0016 


TEST NO.: .... 

OATe: 10 AUGUST 196 2 

COMPUTED by: ?.-Jt 

! p.t u : .123 inch 3 ! n‘ : ' c r u- 3.245 min/in 2 

I STANDPIPE ; LUMP ; mCHES i R I R^ R o 

! AREA I THICKNESS; T n „ cr.TT,Mr. 


TO j op HR I Mlf i j SETTING 


(INCH 2 ) ! (iriCHES) 


.0752 ! 0.05 


CENTERS ! — otT.- 


IN 2 j (INCHES) 


0.025 — 0.251 ±~ .815 T 14.0 


! 8.45 


1.025 2.122 6.89 


5 0.1312 4.242 


.0404 


PAINTED 


0.245 | 0.0016 


1.025 ! 2.122 6.89 


.0756- 0.05 r0.025-T 0.681 4 2.21 “f 11.5 4 


INSIDE SURFACE 
(STILL AIR) 


IR = 5.I75 U =0.1966 



ROOM AIR TEMPERATURE = 24.5 °C 
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APPENDIX A, 


HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 
prori FM- HEAT FLOW THROUGH WALLS test n< 


IME SCALE: I WINr__J HR 


WALL PROFILE 


OUTSIDE SURFACE 
(7.5 MPH WIND) 


I °F r 0.4 INCH 



C“ 

AIR SPACE 



▼ 

FURRING 

i 

0.045 

1.32 

m 

1.2" x 3/4" 

m 

0.101 

1.68 

GYPSUM BOARD 
1/2" 

0.65 

0.45 


INSIDE SURFACE 
(STILL AIR) 


Z R = 5.02 U~.I99 



TEST NOT ; Z 

DATE: 20 _M ARCH 1962 

COMPUTED RY R.H. 

I RTU - 0.0163 5 INCH 3 I °F HR. '^Tll- . 26.10 MIN/IN 2 


STANDPIPE j LUMP j INCHES R p R q R q 

AREA | T H iCKNESS j TO ^ | op HR MIN SETTING 

‘ • CENTERS, i BTu , n 2 (INCHES) 


(INCH* 1 ) ! (INCHES) 


0.0716 1.00 


0.0717 1.00 





0.0728 

1.00 

0.0729 

1.00 

0.0811 

1.25 


5 95 

0.50 I 0.29 7.57 9.85 


1.00 0.08 2.09 11.40 


1.00 ! 0.08 2.09 11.25 


2.09 , 11.30 


0.08 2.09 13.90 


2.09 13.80 


1.00 0.08 2.09 13.70 


1.00 0.08 2.09 11.70 





2.20057.50 


1.970 51.40 6.20 

5.68 
7.85 
5.55 


ROOM AIR TEMPERATURE = 28. 5°C 








































































APPENDIX A. 


HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


ponm 'em; HEAT FLOW THROUGH WALLS 


;E; 20 MARCH 1962 


TIME SCALE 


• | MIN: I HR \ I °F : INCH 


C R 

WALL PROFILE _BTU_ °F H R 

°F BTU 


OUTSIDE SURFACE 
(7.5 MPH WIND) 




COMPUTED BY! fjJl 

I BTU ; 0 01 555 INCH 3 I °V rlR/BTU- 26 .10 MIN/IN 2 


R„ R„ 


) rv o ' 0 

R MIN SETTING 

! Ttu - IN 2 (INCHES) 




0.0720 1.00 


0.0720 I 1.00 


0.0722 


0.0723 


0.0728 | 1.00 


0.0729 1.00 


0.0710 


0.0730 1.00 


0.0730 1.00 


0. 0749 1.00 


0.0311 1.25 


6.05 

0.29 7.57 10.50 


2.09 11.57 



1.00 : 0.08 2.09 11.55 


1.00 ! 0.08 2.09 11.55 


1.00 ! 0.08 2.09 13.98 


2.09 13.88 


1.00 ! 0.06 2.09 13.77 


1.00 ! 0.08 2.09 





1.00 t 0 08 2.09 11.15 


1.00 0.08 2.09 14.02 


1.00 0.08 2.09 14.16 


4.86 

6.95 

1.375 2.408 62.80 1 1| 

15. 15 


0.850 22.40 6 10 


i Rs 4.86 U = .206 


ROOM AIR TEMPERATURE =26.5 W C 
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APPENDIX A 


HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 

pbdri em: HEAT FLOW THROUGH WALLS TEST ni 

date:., 

comput 

rir/.E scale: i win- 1 hr i°F: 0.4 inch irtij : 0.0 468 inch 3 i c e 


C 

BTU 


- TEST NO.: JL 

DAT E 20 MARCH. 1 962 

COMPUTED by: R.H. 

I BTU : 0 0 468 inch 3 i hR 'p.tut ,8.55* min/in 2 


WALL PROFILE 


OUTSIDE SURFACE 
(7.5 MPH WIND) 































































APPENDIX A. 

HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


PROBLEM: HEAT FL OW THROUGH WALLS 


IMF SCALE: I WIN - 1 HR 


I °F 0-4 inch 

-* 1 

! R | 


wa^l ‘ROFile 


OUTSIOE SURFACE 
(7.5 MPH WIND) 

REDWOOD SIDING 
1/2" x 8" LAPPED 

FIR. PLYWOOD 
SHEATING 
( 5/8“ ) 


I JBTU > °F H R 

! °F ! BTU 


0 ! 0.25 


0.726 i 0.85 


1.12 0.78 


-■ TEST NO : !Q- 

OAT e : 20 M ARCH 1962 

COMPUTED RV R-H. 

i RTu= 0-04 80 inch 3 i c f hr/8 7 U - 8 mi n/i n 2 

! STANDPIPE j LUMP j INCHES | R; I R 0 

j AREA | T HICKNESS . TO I o F HR MIN SETTING 

! (INCH 2 ) : (INCHES) | CtNTERS J “btIT - IN 2 (INCHES) 


0.0869' 0.50 


0.1335 i 0.625 


0.0238 1.108 


0.675 5.63 4.70 

9.8 


0.563 | 0.815 6.80 


0.867 : 0.547 4.56 6 10 


4 1.108 i 0.313 2.61 10.5 


0.595 0.94 


0.0240 1.108 


GYPSUM BOARD 
( 1 / 2 ” ) 


INSIDE SURFACE 
(STILL AIR ) 


1.108 * 0.313 2.61 13.1 



ROOM AIR TEMPERATURE = 29.0 • C 













APPENDIX A. 


HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 

ppom fm- HEAT FLOW THR OUG H WALLS t &st ni 

1 da -t e:... 

cow put 

*IMt SCALE: I MIN S .66 6 hr I INCH !BTU-__^P267 inch 3 t °F 


AALL PROFILE | 


CONCRETE 


i BTU 
! 


I '’F - INCH 

1 1 

! R ! 

i _°f HIT | 

I BTU 


T&ST no - PRO BLEM I 

DA ’e: 12 JUNE 1962 

COMPUTED by: B. L_. 

l°F HR P.TU- .01775 min/in 3 

“ ! r„ I rT ! rT ~1 


i STANDPIPE j LUMP |NCHES I R R R 

; AREA ITHICKNESS, Tn j „ cr TTi«r 

i I IO j Of- gg MIN SETTING 

! (INCH 2 ) (INCHES) 1 CENTERS | "gfQ ~ |N 2 (INCHE S) 

! i 

I .0400 i .80 L - 4 032 1 DU: 

! — j — j .804 ! .0643 3.62 13.05 

i .0404; . .808 | — - — ■ — f— 

1 1 i .814 ; .0651 3.66 8.45 

| .0410 j .82 i 1 

| 1 1 .828 ; .0663 3.72 8.55 

! .0418 ! .836.1- 

.838 : .067 3.76 12.07 

.848 i .0678 3.81 13 04 

1.103 .0883 4.96 12.77 

— i — ! — u — 

1.371 I .1098 6.17 4.6 

1.406 | .1124 j 6.32 16.3 

i 

! f ■ 

■ i 

1.43 ! .1144 i 6.43 II. I 


.0420 ! .84 

i 

J 

0428 ! .856 


j 

.0675 

— 

1.35 

! 

.0696 

— 

i 

1.392 ! 

j 


.0710 1 .42 


.0720! 1.44 


.0720 1 1.44 


1.44 

! .1151 

6.48 

.72 

1 .0576 

3.24 


INSULATION 
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HYDRAULIC ANALOGUE COMPUTER SCALING SHEET 


ppnRi em: heat flow through walls 

TIME SCALE: I MIN: ! HR I °F : 5/8 _ inch 


-J 

wall profile 

C 

BTU 

°F 

1 

R 

°F HR 
BTU 




MATERIAL "X" 

11.25 

0.96 

INSULATION 




test mi v PROBLEM 2 
DAT E ; 21 JULY 1962 


COMPUTEO by: R H. 

IBTU-- INCH 3 l°r HR/BTUr 26.05 ,_MIN/!N 2 


standpipe 

AREA 
(INCH 2 ) 

LUMP 

THICKNESS 

(INCHES) 

INCHES 

TO 

CENTERS 

R p 

°F HR 
BTU 

R o 

MIN 

IN 2 

R c 

setting 

(INCHES) 


1 

i 





0.0720 



0.50 j 0.08 

\ 

2.08 

11.05 

1.00 

1 

1 

i nn 

ft 1C 

4.16 

7.73 

0.0722 


I ' 

i 

1 .00 

. 

t 

i 

i 

1.00 | 0.16 

1 

1 

i 

4.16 

7.05 

0.0720 

1.00 

— 1 

i 

j 

■ AA ! A 1C 

4.16 

7.17 

0.0723 

1.00 

1 .ww 

■ 

V/. IO 

1.00 

0.16 

4.16 

12.5 

0.0717 

1.00 

1.00 

0.16 

4.16 

11.2 

0.0716 

1.00 

0.50 

0.08 

2.08 









ROOM AIR TEMPERATURE" 28.0°C 



APPENDIX B. 


DETERMINING PARAMETERS IN THE THERMAE NETWORK 

It has been shown that the replacement of a field problem by a network approxi- 
mation, such as the thermal network used in this investigation, leads to an inevitable 
error in the final result. This error will be 5% or less if the thickness of the section 
is chosen such that it is equal to or less than one-eighth of a thermal wavelength 4 : 



A x < L = 0. 44 ft (Bl) 

where L = one-eighth thermal wavelength (X./8), ft. 

/ k \ , \ 
a = thermal diffusivity (— — ■), ft 2 /hr. 

P P 7 

f = frequency, cycles/hr. 

In Table II the maximum mesh sizes for an error of 5% or less are recorded for the 
construction materials used in the ten walls, taking f = 1/24 cycles/hr. It is observed 
that all of the mesh sizes used in the ten walls tested were considerably smaller than 
the values listed in this table. 

From eq Bl above, it is seen that, if the frequency is increased, the mesh size 
necessary for the same degree of accuracy is decreased. Since the input curve is 
composed of a constant plus a diurnal sine wave and several higher frequencies signifi- 
cant to the fourth harmonic, sections about half the size given in Table II are needed to 
hold the 5% tolerance. Each section should contain only one material, where possible, 
for this estimate of the error to be accurate. 


APPENDIX C. 

METHOD USED IN DETERMINING THE TIME WHEN THE INPUT 
WAS CYCLING PERIODICALLY IN THE STEADY STATE 

In any ladder network, like the thermal network pictured in Figure 2, the input 
must cycle for a certain period of time before the system is in the steady state and 
the transients have died out. A glance at the thermal network reveals that it is com- 
posed of several loops with several different time constants. The time necessary for 
the system to reach equilibrium is dependent upon the loop wherein the time constant 
is la'rgest: / 


Time constant - RC = hr. 
where 0=0! + C 2 . 


For every time constant elapsed, 63% of the transients still in the loop will have dis- 
appeared, or 99. 5% will have disappeared after the input has been cycling for five 
time constants. This time was converted to computer time — the computer was allowed 
to run for this length of time plus half as much again before data were taken. This 
time rarely exceeded three cycles and was usually much less. 



APPENDIX D. 


CHECK ON HYDRAULIC ANALOG COMPUTER ACCURACY 

The hydraulic analog computer was tested in the first phase of the investigation 
by application of a sine wave and step input to the outside surfaces of a homogeneous 
slab of width 2 i (see test problems 1 2). The results from the computer were then 

checked against the analytical solutions 5 (Fig. Dl, D2). In Figure D1 (the step input), 
the difference between the analog and analytical solutions is less than 1% for 

0' 1 ' > 0. 10. When 0' 1 ' = 0. 005,. the computer value is about 10% lower than the 

analytical value. The latter deviation can be satisfactorily explained by realizing that 
the flow in the hydraulic resistors is turbulent whenever the head is- much larger than 
3 in. (The head for the first few seconds of the step problem was almost 5 in. ) 
Turbulence increases the hydraulic resistance, which results in a decrease of current 
and potential. This problem was not encountered in testing any of the ten wall sections. 

The sine wave problem (Fig. D2) shows reasonable agreement with the analytical 
solution, although the error in the time lag is of the order of 4 or 5%. This error was 
found to be due largely to the input being fed by the kymograph at a speed slower (3-4%) 
than expected. This error in speed had the effect of increasing the effective scale of 
resistance. The net result, is to lower the amplitude and change the time lag. 

The accuracy of the computer was therefore considered reasonable and the sources 
of error in these two test problems were corrected as much as possible before studying 
the actual wall problems. 

DEFINITION OF TERMS 

temperature 
period of the function 
2 tt/P 

phase angle sine wave input temperature 
phase angle of v within the solid 
thermal conductivity of solid 
specific heat of solid 
density of solid 
thermal diffusivity of solid 
time 

Test problem 1. Step function. 

0 < x < ZH 

v = v (x, t) 

Governing equation: 

k _ 9 Z v _ _9v 
pCp ’ ax* " 8t 

0 l Zl 

X » 



V — 

p - 
u — 
€ — 
<t> — 
k - 
Cp - 

P - 
a, a — 
t — 


APPENDIX D 


LEGEND 

ANALYTICAL 

o- o COMPUTER 


9 *= /. 50 


\\\X V v 
\\ X \X X . 

v X \ \\\ 


' Y\\ \\ 


'\ \ \ \\ X X x 


A \ \ \\ X x X \ X x 

A\ V \ V ' 

A \ V \x \ x x0 '? 

\ \ \ \V xo X ' 

X \ \ \\Vo 

\ \ Vp ^ 

\ \ %\\ 

\p Xx "X 


> \ 

XX ' 


X 

. x<p 
^ . o . x - 


Figure Dl. Temperature variations in a uniform slab. 


LEGEND 


XX 


\\\ 


PEAK TO PEAK AMPLITUDE OF SURFACE V 


AVERAGE V IN STEAOY STATE 


ANALYTICAL 


-O COMPUTER 


W\ 


\\\ 


Figure DZ. Temperature history at two depths in a 
homogeneous solid with harmonic surface temperature. 
















APPENDIX D. 


D3 


Boundary conditions: 

v (21 , t) = v (o, t) = V t > 0 

v (x, t) = 0 t < 0 

v (x, 0) = 0 for 0 < x < 2f 



at x - i 


Define : 


Y _ ^ , "T ' V . 

x ~ t ’ 1 v ’ 


a 


k 

pCp ’ 


e* 


at 

f 2 


9 2 T _ 3T 

ax 2 a e" 


Solution: 



V = A sin (cot + e + 4») 


+ 4 




n = 0 


( - if 1 ( Zn + 1) [ 4l 2 cj cos e - k (2n + l) 2 ir 2 sin e] 
1 6 St 4 u> 2 + k 2 it 4 (2n + I) 4 


-k(2n + l) 2 it 2 t/4f 2 (2n + 1 ) tt x 

e _ cos u 


t Car slaw & Jaeger, Conduction of heat in solids, Second Edition, 
Clarendon Press, 1959> p. 100. 

** Ibid. , p. 105 


wave. 

-I < x < t 
x > |i | 
t > 0 


(steady state) 


(transient)*^ 


D4 . APPENDIX D. 

Assume steady state reached: 

V = A sin (cot + e + ©) 

_ j cos h y X ( 1 + i) | _ f cos h 2 yx + cos 2~y x j 2 
cosh \t (1+i) ,cosh 2 yi + cos i / 


4, = ARC \x(l til l 

Y ! co sh yi ( 1 + 1 ) l 


y = ( to/ 2 a ) 2 
a = k/ pCp 
P = 2 ir/to 

Calculations 

x — 0. 292 ft 


£ = 0. 500 ft 


to = 0. 2 62 /hr 


k = 0. 521 Btu/hr -ft- ° F 
a = k/pCp = 0.0231 ft 2 /hr 
pCp = 22. 50 Btu/ft 3 -°F 

l 

Y = (to/ 2 a) 2 = 2. 38/ft 

e = 90° v = A cos (cot + <j>) 

/ 1 
A _ / cos h 1. 39 + cos 78. 6 \ z _ n 

A “ Uosh 2.38 + cos 136. 3/ “ °‘ 696 

0 =ARG .Jc°sMx.(l+i)1 
(cos h yl (1 + i )J 

cosh Y x(l +i) = He VX(1+i) +e' vx(1 + i) ] 

= | [e ' x (cos yx + i sin yx) + e ^(cosyx - i sin yx)] 
= | [cos yx (2 cosh Y x ) + i sin yx. (2 sin h y x )] 

9^ = tan -1 |tan y x tanh y x ^ - 26° 37' 

0q = tan -1 "I tan y£ tanh y£ "j- = 64°20' 

4> = °n - °D = - 37 ° 43, 


x = 0. 125 ft 


f cosh 0. 595 + cos 34. 1 
'.cosh 2. 38 + cos 13 6. 3 


= 0. 650 


0 N =5"5' 

0 D = 64° 20' 

4> = - 59° 15'-. 


(other parameters same) 


APPENDIX E. 


CHECK ON HEAT FLOW RESULTS 

As an approximate check on the results obtained for the heat flow through the 
ten wails, the thermal circuit was divided into two sections (Table III) and a sine 
wave with a frequency of 1/24 cycles/hr was used as the input: 



f( t) - A 4 B sin( 2 tt ft) P. F. = Peak ioad factor 


A = 20° F 


A/R + B/Z 
A/R ' 


B = 8 ° F 


0 = Time lag, hrs. 


Values for Z and 0 for the ten walls were obtained by a method outlined by Guillemin 6 . 
The resulting time lags and peak load factors provided by this analytical approximation 
were recorded in Tables IV and V, respectively. 


APPENDIX F. 


PHOTOGRAPHIC TECHNIQUE FOR RECORDING DATA 

Most of the wall sections tested in this investigation were divided into 8 to 14 
sections. It was decided that to record the height of the fluid manually in ten or more 
standpipes every minute (analog equivalent of one hour in the thermal system) would 
be both impracticable and inaccurate. Consequently, the authors chose a method of 
recording data which consists of photographing these heights. 

A 35 millimeter camera was set up in front of the computer, at a distance just 
great enough to encompass all of the standpipes. A string of five flood lights was 
placed about 2 ft behind the standpipes. A large sheet of onion skin paper was placed 
directly in back of the standpipes to diffuse the light. To make the fluid show up 
against the light background, ordinary food coloring was added to the distilled, de-aired 
water. (A check was made to see whether the food coloring affected the viscosity of 
the water. The change in viscosity was negligible.) The standpipes were coated with 
Drifilm* to eliminate the meniscus and 0. 01% mercuric chloride solution was added to 
the water to prevent the formation of bacteria. In order to read the heights accurately 
from the film, two horizontal refer ence lines placed 1 5 in. apart were set up directly 
in front of the standpipes. 

Pictures were taken at one -minute intervals over the duration of one cycle. The 
films were developed in the laboratory. To read the heights, the film was projected 
on a large screen which had already been calibrated to tenths of an inch. 

In the two test problems, the error produced by this method of recording data 
was found to be less than 1%. The photographic method was subsequently used in the 
investigations of the wall sections. 


* Drifilm 0087. General Electric Organosilicon Products 


